Day-night effect in solar neutrino oscillations with three flavors 
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We investigate the effects of a nonzero leptonic mixing angle On on tlie solar neutrino day-night 
asymmetry. Using a constant matter density profile for the Earth and well-motivated approxima- 
tions, we derive analytical expressions for the Ve, survival probabilities for solar neutrinos arriving 
directly at the detector and for solar neutrinos which have passed through the Earth. Furthermore, 
we numerically study the effects of a non-zero ^13 on the day-night asymmetry at detectors and find 
that they are small. Finally, we show that if the uncertainties in the parameters 612 and Am'^ as 
well as the uncertainty in the day-night asymmetry itself were much smaller than they are today, 
this effect could, in principle, be used to determine On- 

PACS numbers: 14.60.Pq, 13.15.+g, 26.65. +t 



I. INTRODUCTION 

Neutrino oscillation physics has entered the era of pre- 
cision measurements with the results from the Super- 
Kamiokande Q, SNO 0, S S H, and KamLAND Q 
experiments. Especially, impressiver results have re- 
cently come from measurements of solar neutrinos (see 
Refs. 0, 0, 0) and the solar neutrino problem has suc- 
cessfully been solved in terms of solar neutrino oscilla- 
tions. 

The solar neutrino day-night effect, which measures 
the relative difference of the electron neutrinos coming 
from the Sun at nighttime and daytime, is so far the best 
long baseline experiment that can measure the matter ef- 
fects on the neutrinos, the so-called Mikheyev-Smirnov- 
Wolfenstein (MSW) effect In all accelerator long 

baseline experiments, the neutrinos cannot be made to 
travel through vacuum. The atmospheric neutrino exper- 
iments, on the other hand, use different baseline lengths 
for neutrinos traversing the Earth and those that pass 
through vacuum. With the advent of the precision era 
in neutrino oscillation physics, we can gradually hope to 
obtain better measurements of the day-night effect. Re- 
cently, both Super-Kamiokande and SNO have presented 
new measurements 0, of this effect that have errors 
approaching a few standard deviations in significance. 

In this paper we analyze the day-night effect in the 
three neutrino flavor case. Earlier analyses of this ef- 
fect, with a few exceptions, have been performed for 
the two neutrino flavor case. Furthermore, the present 
data also permit a new treatment of the effect due 
to the particular values of leptonic mixing angles and 
neutrino mass squared differences obtained from other 
experiments. There are six parameters that describe 
the neutrinos in the minimal extension of the standard 



model: three leptonic mixing angles 9i2, O13, and 6*23, one 
CP-phase S, and two neutrino mass squared differences 
AM^ = 7713 — ml and Am^ = — . The solar neu- 
trino day-night effect is mainly sensitive to the angles 9i2 
and 013, and the mass squared difference Am^. Our goal 
is to obtain a relatively simple analytic expression for the 
day-night asymmetry that reproduces the main features 
of the situation. It turns out that one can come a long 
way towards this goal. 



Earlier treatments of the day-night effect can be found 
in Rcfs. [i mill mill El lli. Our three flavor treat- 
ment is consistent with the modifications presented by 
de Holanda and Smirnov |l2j | as well as Bandyopadhyay 
et al. 13. 
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This paper is organized as follows. In Sec. ^ we in- 
vestigate the electron neutrino survival probability with 
n flavors for solar neutrinos arriving at the Earth and 
for solar neutrinos going through the Earth. Next, in 
Sec. mil we study the case of three neutrino flavors, in- 
cluding production and propagation in the Sun as well as 
propagation in the Earth. At the end of this section, we 
present the analytical expression for the day-night asym- 
metry. Then, in Sec. II VI we discuss the day- night effect at 
detectors. Especially, we calculate the elastic scattering 
day-night asymmetry at the Super-Kamiokande experi- 
ment and the charged-current day-night asymmetry at 
the SNO experiment. Furthermore, we discuss the possi- 
bility of determining the leptonic mixing angle ^13 using 
the day-night asymmetry. In Sec. we present our sum- 
mary as well as our conclusions. Finally, in the Appendix 
we shortly review for completeness the day-night asym- 
metry in the case of two neutrino flavors. 
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II. THE n FLAVOR SOLAR NEUTRINO 
SURVIVAL PROBABILITY 

Assuming an incoherent neutrino flux [ToL ITll | , the Vf. 
survival probability for solar neutrinos is 

n n 
i=l i=l 

where n is the number of neutrino flavors and ki is the 
fraction of the mass eigenstate \ui) in the flux of solar 
neutrinos. From unitarity it follows that 



(2) 



1/n for all i, we 



In the case of even mixing, i.e., k 
obtain Ps = 1/n. 

For neutrinos reaching the Earth during daytime (at 
the detector site), Ps is the i^e survival probability at the 
detector. However, during nighttime this survival prob- 
ability may be altered by the influence of the effective 
Earth matter density potential. Thus, in this case, the 
survival probability becomes 



P.' 



SE 



Y^k. 

4=1 



(3) 



where \vi) = \vi{L)) and L is the length of the neutrino 
path through the Earth. Here, the components of \vi{t)) 
satisiy the Schrodinger equation 



dt 



H,n\Mt))r 



(4) 



and where 



with the initial condition |^'i(0)) = 
denotes the mass eigenstate basis. 

The Hamiltonian H is given by (assuming k sterile 
neutrino flavors) 



k times 



i/„, ~ — + C/Tdiag(ycc,0,. 
2p 



.,0,-Vmc,---,-Vnc)U, 
(5) 

where M — diag(mi, m2, . . . , rUn), the effective charged- 
current Earth matter density potential is Vcc = 
V^Gf^c, and the effective neutral-current Eath mat- 
ter density potential is Vnc = where Gp 
is the Fermi coupling constant and where and iV„ are 
the electron and nucleon number densities, respectively. 
The number densities are functions of t depending on 
the Earth matter density profile, which is normally given 
by the Preliminary Reference Earth Model (PREM) [l6l |. 
The term | {v^lvi) P is interpreted as the probability of a 
neutrino reaching the Earth in the mass eigenstate {vi) to 
be detected as an electron neutrino after traversing the 
distance L in the Earth. For notational convenience we 
denote 



P. 



(6) 



Clearly, P,,{L = 0) = \ 
from unitarity it follows that 



|[/eiP. Furthermore, 



(7) 



i=l 



Again, in the case of even mixing, we obtain Pse = 
and the Ve survival probability is unaffected by the pas- 
sage through the Earth. 



III. THE CASE OF THREE NEUTRINO 
FLAVORS 

Until now most analyses of the day-night effect have 
been done in the framework of two neutrino flavors. How- 
ever, we know that there are (at least) three neutrino fla- 
vors. The reason for using two flavor analyses has been 
that the leptonic mixing angle ^13 is known to be small 
|l7j . leading to an approximate two neutrino case. One 
of the main goals of this paper is to find the effects on the 
day-night asymmetry induced by using a non-zero mix- 
ing angle ^13. In what follows, we assume that there are 
three active neutrino flavors and no sterile neutrinos. 

We will use the standard parametrization of the 3x3 
leptonic mixing matrix |l8j| 



U 



where s,,- = sin( 



C13C12 C13S12 si3e 

* * * 

* * * 



(8) 



and 



cos 0ij, 0ij are leptonic 



mixing angles, and the elements denoted by * do not 
affect the neutrino oscillation probabilities, which we are 
calculating in this paper. 



A. Production and propagation in the Sun 

In the three flavor framework, there are a number of 
issues of the neutrino production and propagation in the 
Sun, which are not present in the two flavor framework. 
First of all, the three energy levels of neutrino matter 
eigenstates in general allow two MSW resonances. Fur- 
thermore, the matter dependence of the mixing parame- 
ters are far from as simple as in the two flavor case. The 
result of this is that wc have to make certain approxima- 
tions. 

Repeating the approach made in the two flavor case 
(see the Appendix), we obtain the following expression 
for h: 



dr/(r 



3 

E 



(9) 



where U is the mixing matrix in matter, P^^ is the proba- 
bility of a neutrino created in the matter eigenstate {i^aij) 
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to exit the Sun in the mass eigenstate \ vi), and /(r) is the 
normaUzed spatial production distribution in the Sun. 

The second resonance in the three flavor case occurs 
at Vcc — cos{29i3)AM-^ /{2E), assuming that the reso- 
nances are fairly separated. The maximal electron num- 
ber density in the Sun, according to the standard solar 
model (SSM) [l9l|, is about iVe,,„ax ^ 102 TV^/cm^, yield- 
ing a maximal effective potential Vccmax — 7.8 x 10"^^ 
MeV. Assuming the large mass squared difference AM^ 
to be of the order of the atmospheric mass squared differ- 
ence (jAm^t^l ~ 2 X 10~^ eV^ US): the neutrino energy 
E to be of the order of 10 MeV, and 6*13 to be small, we 
find |AMV2i;|cos(26'i3) ~ IQ-^^ ^gy. Thus, the solar 
neutrinos never pass through the second resonance, inde- 
pendent of the sign of the large mass squared difference. 

Since the neutrinos never pass through the second res- 
onance, it is a good approximation to assume that the 
matter eigenstate {vm.s) evolves adiabatically, and thus, 
we have 

P3k^Pk3=S3k- (10) 

Unitarity then implies that 

P12 ^ P2I = -Pjump- (11) 

Furthermore, if we assume that Vcc ^ Am'^/{2E) <C 
AM^/(2i?), the neutrino evolution is well approximated 
by the energy eigenstate \vm,3) evolving as the mass 
eigenstate I1/3) and the remaining neutrino states oscil- 
lating according to the two flavor case with the effec- 
tive potential VcS — cf^Vcc- This does not change the 
probability Pjump as calculated with a linear approxima- 
tion of the potential in the two flavor case. This means 
that we may use the same expression as that obtained 
in the two flavor case, see the Appendix, even if the res- 
onance point, where |A^e/-^e| is to be evaluated, does 
change. However, in the Sun, is approximately expo- 
nentially decaying with the radius of the Sun, leading to 
\Ne/Ne\ being approximately constant, and thus, inde- 
pendent of the point of evaluation. For the large mixing 
angle (LMA) region, the probability -Pjump of a transition 
from \vm,i) to \vm,2), or vice versa, is negligibly small 
(Pjump < 10^^''"°). However, we keep it in our formulas 
for completeness. 

To make one further approximation, as long as the 
Sun's effective potential is much less than the large mass 
squared difference AM^, the mixing angle ^13 ~ ^13 giv- 
ing 

fc3~ / dr/(r)sin2 0i3 =sin2 0i3. (12) 

JQ 

A general paramctrization for ki and ^2 is then given 
by 

ki^c^s ^ , k2^c^3 . (13) 

In the above approximation, the oscillations between the 
matter eigenstates |i^a/.i) and |i^m,2) are well approxi- 
mated by a two flavor oscillation, using the small mass 



difference squared Am?, the mixing angle 612, and the 
effective potential Vcs = c\^Vcc- Thus, we obtain 

D3,= / dr/(r)cos[2^i2(r)](l-2Pjump), (14) 
Jo 

where cos[26'i2(r')] is calculated in the same way as in the 
two flavor case using the effective potential. For reason- 
able values of the neutrino oscillation parameters, this 
turns out to be an excellent approximation. 

Inserting the above approximation into Eq. (^) with 
n = 3, we obtain 

p 4,4 I + -D31. C0S(26>12) 

Ps = Si3 + Ci3 . (15) 

When 6*13 0, we have D31, D^v, and thus, we recover 
the two flavor survival probability in this limit. 



B. Propagation in the Earth 

As in the case of propagation in the Sun, Vcc ^ 
Am?/2E < AM'^/2E, \vm,3) ^ Ws), and the remaining 
two neutrino eigenstates evolve according to the two fla- 
vor case with an effective potential of VcS — cf^Vcc- For 
the MSW solutions of the solar neutrino problem along 
with the assumption that AM^ is of the same order of 
magnitude as the atmospheric mass squared difference, 
this condition is well fulfilled for solar neutrinos propa- 
gating through the Earth. As a direct result, we obtain 
the probability P^e as 

P3, ~|(z.,|z.3)P=|C/e3|' = 4- (16) 

It also follows that 

P2e = 4^ sin2(20i2) sm\aL) + cj^sj^, (17) 
where 

a^^^K^- 2CI3VEK cos{2ei2)+ct3Vi, (18) 

Ve is the electron neutrino potential in the Earth, and 
K = Am} I2E. We observe that when i = or Ve = 0, 
Pze = (^3^X2 — I^e2p just as cxpectcd. 

C. The final expression for Pn-d 

Now, we insert the analytical expressions obtained in 
the previous two sections into Eq. (^J with n — 3 and 
subtract P5 from this in order to obtain an expression for 
Pn-d = PsB ~ Ps in the three flavor framework. After 
some simplifications, we find 

Pn-d = -cl^Ds,^ sin2(20i2) sin2(aL). (19) 
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For the MSW solutions of the solar neutrino problem, 
K > c^gVe, and thus, a ~ K/2. This yields 

P„_d ^ -24^3. 1^ sin2(20i2) sin^ (^l) . (20) 

Apparently, the effect of using three flavors instead of 
two is, up to the approximations made, a multiplication 
by as well as a correction in changing D21, to D^^. 
When 0, we have Z?3y £'21/, and we regain the 

two flavor expression in this limit (see the Appendix) . An 
important observation is that the regenerative term, for 
Ve <C K, is linearly dependent on Ve- Thus, the choice 
of which value of Ve to use is crucial for the quantitative 
result. As is argued in the appendix, the potential to 
use is the potential corresponding to the electron number 
density of the Earth's crust. However, the qualitative 
behavior of the effect of a non-zero 6*13 is not greatly 
affected. 
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FIG. 1: The zenith angle exposure function Y{a) for SK and 
SNO as a function of the zenith angle a. The data are re- 
trieved from Ref. '21]. 



IV. THE DAY-NIGHT EFFECT AT 
DETECTORS 



From the calculations made in the previous parts of 
this paper, we obtain the day-night asymmetry of the 
electron neutrino flux at the neutrino energy E as 



AtuiE) = 2 



C^e.NjE) - <l>e,D{E) 
(l)e,N{E) + (t>e^D{E) 
Pn-d{E) 



Ps{E) 



-aiE) 



(21) 



However, this is not the event rate asymmetry measured 
at detectors. We will assume a water-Cherenkov detector 
in which neutrinos are detected by one of the following 
reactions |31J: 



P + P + e" 



(22) 
(23) 



where x = e,fj,, r, which are referred to as elastic scat- 
tering (ES) and charged-current (CC), respectively. The 
CC reaction can only occur for x = e, since inserting 
a; ^ e in Eq. I|23|) would violate the lepton numbers Le 
and Lx- We assume that the scattered electron energy 
T' is measured and that the cross sections clcr^^^ / dT' and 
da^^/dT' are equal. 

If we denote the zenith angle, i.e., the angle between 
zenith and the Sun at the detector, by a, then the event 
rate of measured electrons with energy T in the detector 
is proportional to 



i?(a,T) = / dE(j){E) 

Jq Jo 

xF{T,T''^'' 



dT' 



dT' 



(24) 



where (ji{E) is the total solar neutrino flux, T' is the true 
electron energy, and da^^^^.^JdT' is given by 



da,. 



dT' 



PSE 



d(T„ 



dT' 



(I-Pse)' 



da,. 



dT' 



(25) 



Here, we have used the assumption dcTj^^/dT' — 
d(Ti/^/dT', since neutrinos not found in the state lue) are 
assumed to be in the state or in the state Ifr)- The 
energy resolution of the detector is introduced through 
F{T,T'), which is given by 



F{T, T') 



exp 



(T - T'f 
2A|r" 



(26) 



where A^' is the energy resolution at the electron energy 

r. 

The night and day rates N and D at the measured 
electron energy T are given by 



D{T) = / daR{a,T)Y{a), 
Jo 



N{T) 



Tr/2 



daR{a,T)Y{a), 



(27) 
(28) 



respectively. Here, Y(a) is the zenith angle exposure 
function, which gives the distribution of exposure time 
for the different zenith angles. The exposure function 
is clearly symmetric around a — 7r/2 and is plotted in 
Fig. ^for both Super-Kamiokande (SK) and SNO. From 
the night and day rates at a specific electron energy, we 
define the day-night asymmetry at energy T as 



An-d{T) = 2 



N{T) - D{T) 
N{T) + D{T)' 



(29) 



The final day-night asymmetry is given by integrating the 
day and night rates over all energies above the detector 
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threshold energy Tth, i.e., 



A 



J^JT[NiT)~DiT)] ^N-D 



n—d 



= 2- 



/^^ dr[iV(T) + D{T)] N + D 



(30) 



The threshold energy Tth is 5 MeV for both SK and SNO. 

For computational reasons, we will start by performing 
the integral over the zenith angle a. For the daytime flux 
D, PsE — Ps, which is independent of a. As a result, 
the only a dependence is in Y{a) and the zenith angle 
integral only contributes with a factor one-half [if the 
normalization of Y is such that AaY{a.) = 1]. In order 
to be able to use the results we have obtained for P„_d, 
we need to compute the difference between the night and 
day fluxes, which is given by 



N{T) - D{T) = / daY{a)R 

Jtt/2 



(31) 



where the quantity i?" — R{a, T) — i?(7r — a, T) is on 
the form 

n-d _ I ATP At TP \ I fYT' 



Jo Jo 

xi?(T,r')-^ 



(32) 



and 



— = Pn-dia, E,)[-^- ) . (33) 



dT' 



V dT' dT' 



Note that the a dependence in Pn-d enters through the 
length traveled by the neutrinos in the Earth and that 
the argument aL of the second sin^ factor in Eq. H19|) 
oscillates very fast and performs an effective averaging of 
Pn-d in the zenith angle integral, i.e., replacing sin^(aL) 
by 1/2. After this averaging, the only zenith angle depen- 
dence left is that of Y{a) and the zenith angle integral 
only gives us a factor of one-half as in the case of the day 
rate D. 



A. Elastic scattering detection 

Neutrinos are detected through ES at both SK and 
SNO. The ES cross sections in the laboratory frame are 
given by Ref. [3- For kinematical reasons, the maximal 
kinetic energy of the scattered electron in the laboratory 
frame is given by 



T' 



1 



(34) 



2E^ 



The integrals that remain cannot be calculated ana- 
lytically. Hence, we use numerical methods to evaluate 
these integrals. However, computing all integrals by nu- 
merical methods demands a lot of computer time, and 
thus, we make one further approximation, that all solar 




FIG. 2: The day-night asymmetry at SK for different val- 
ues of T, Am^, and Orz as a function of 9\3 relative to the 
corresponding value for ^13 = 0. 



neutrinos are produced where the solar effective po- 
tential is Vcc ^ 7.07 X 10"!*^ MeV, which is the effective 
potential at the radius where most solar neutrinos are 
produced. For reasonable values of the fundamental neu- 
trino parameters, the error made in this approximation 
is small. 

For the energy resolution of SK, we use ^3 



At' = 1.6 MeVv/r'/(10 MeV), 



(35) 



and for the electron number density in the Earth, we 
use Nf, = IANa/chi^, where Na is the Avogadro con- 
stant, which roughly corresponds to 2.8 g/cvn^ (using 
Z/A ~ 0.5, where Z is the number of protons and A 
the number of nucleons for the mantle of the Earth). 
The electron number density used corresponds to the 
density in the Earth's crust. The motivation for using 
this density rather than a mean density can be found in 
the Appendix. Note that the regenerative term Pn-d in 
Eq. (|20|) . and thus, the day- night asymmetry, is linearly 
dependent on the matter potential Ve- It follows that the 
electron number density used has a great impact on the 
final results. If we had used the average mantle matter 
density of about 5 g/cm'^, then the resulting asymmetry 
would increase by almost a factor of two. 

The above values give us the numerical results pre- 
sented in Fig. 13 As can be seen from this figure, the rela- 
tive effect of a non-zero ^13 is increasing if the small mass 
squared difference Am^ increases or if the measured elec- 
tron energy or the leptonic mixing angle 9i2 decreases. 
The effect of changing O12 is also clearly larger for smaller 
electron energy T and larger small mass squared differ- 
ence Ato^. In Fig.Othe isocontours of constant day-night 
asymmetry in the SK detector with ^13 equal to 0, 9.2° 
and 12° are shown for a parameter space covering the 
LMA solution of the solar neutrino problem. The val- 
ues used for ^13 correspond to no mixing as well as the 
CHOOZ upper bound for AM^ equal to 2.5 x lO'^ eV^ 
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FIG. 3: Isocontours in the parameter space for the 

ES day-night asymmetry for three different values of Q\z. The 
values of A„-d for the different isocontours are shown in the 
figure. The shaded regions correspond to the allowed regions 
of the parameter space for different confidence levels and the 
circle corresponds to the best-fit point according to Ref. [20|. 



FIG. 4: The CC day-night asymmetry at SNO for difiterent 
values of T, Am^, and 612 as a function of 613 relative to the 
corresponding value for ^13 = 0. 



and 2.0 X lO^'^ eV^, respectively can be seen from 

this figure, the variation in the isocontours are small com- 
pared to the size of the LMA solution and to the current 
uncertainty in the day-night asymmetry 0,0]. However, 
if the values of the parameters 9i2, Am^, and A„_rf were 
known with a larger accuracy, then the change due to 
non-zero ^13 could, in principle, be used to determine 
the "reactor" mixing angle ^13 as an alternative to l ong 
baseline experiments such as neutrino factories l24| 
and super-beams 23^ 25^, as well as future reactor 
experiments and matter effects for supernova neutri- 
nos . The day-night asymmetry for the best-fit value 
of Ref. |23| is An-d — 3.0 %, which is larger than the 
theoretical value quoted by the SK experiment, but still 
clearly within one standard deviation of the experimental 
best-fit value A,,^d = (1-8 ± l-GtH) % 0- 

B. Charged-current detection 

Only SNO uses heavy water, and thus, SNO is the only 
experimental facility detecting solar neutrinos through 
the CC reaction (^5)1 . Since the electron mass me is 
much smaller than the proton mass rup (rUe ~ 511 keV, 
nip ~ 938 MeV) , most of the kinetic energy in the center- 
of-mass frame, which is well approximated by the labo- 
ratory frame, since the deuteron mass by far exceeds the 
neutrino momentum, after the CC reaction will be car- 
ried away by the electron. This energy is given by 

T' = ^ + AS„,a.s, (36) 

where Ai^mass = md — 2mp — me ~ —1.95 MeV. Thus, we 
approximate the differential cross section dcr^^ /dT" by 

^ ^a,J{T' -E +1.95 MeY). (37) 



In the above expression, we use the numerical results 
given in Ref. p3 | and perform linear interpolation to cal- 
culate the total cross section a,^^ as a function of the 
neutrino energy E. For a; 7^ e, the reaction in Eq. H23|) 
is forbidden, since it violates the lepton numbers Lg and 
Lx- Thus, for x ^ e, we have da,^^/dT' ~ 0. 

The energy resolution at SNO is given by 0, 



At' = -0.0684 MeV 

-h0.331 MeV ^(T'/MeV) 

+0.0425 MeV (T'/MeV). (38) 

This gives the results presented in Fig. ^ This figure 
shows the same main features as Fig. |21 However, the 
effects of different T and Am^ are larger in the CC case. 

In Fig [SI we have plotted isocontours for the CC day- 
night asymmetry for 6*13 equal to 0, 9.2° and 12° in order 
to observe the effect of a non-zero 6*13 for the day-night 
asymmetry isocontours in the region of the LMA solu- 
tion. Just as in the case of ES, the isocontours do not 
change dramatically and the change is small compared 
to the uncertainty in the day-night asymmetry. It is also 
apparent that the day-night asymmetry is smaller for the 
ES detection than for the CC detection. This is to be ex- 
pected, since the ES detection is sensitive to the fluxes of 
Vfj, and Vr as well as to the Ve flux, while the CC detection 
is only sensitive to the Ve flux. The day-night asymme- 
try for the best-fit values of Ref. ^ is An-d — 4.7%, 
which corresponds rather well to the values presented in 
Refs. ^3i^3- The latest experimental value for the day- 
night asymmetry at SNO is An-d = (7.0 ± 4.9^1-2) % 
4]. 
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FIG. 5: Isocontours for the CC day-night asymmetry in the 
6i2-Am^ parameter space for three different values of ^13. 
The values of An-d for the different isocontours are shown in 
the figure. The shaded regions and the circle are the same as 
in Fig. El 



values presented give a notion about their magnitudes. 
If we use a more optimistic estimate instead of Eq. (|39|l , 
then S9i3 reduces to about 8° . Unfortunately, this is still 
a quite large uncertainty and in order to reach an uncer- 
tainty of about 1°, we need measurements of An-d and 
Am^ with uncertainties that are about 100 times smaller 
than the uncertainties of today and a measurement of 0i2 
with an uncertainty that is about 10 times smaller than 
today. 

We would like to stress that the above can only be 
considered as an indication of the uncertainties needed 
to determine 613 and serves mainly to present the idea 
of measuring ^13 with the solar day-night asymmetry. In 
a realistic scenario, uncertainties and fluctuations in the 
Earth matter density must be treated in a more rigorous 
way than using a constant Earth potential; this would 
probably require a full numerical simulation with three 
neutrino flavors. 



V. SUMMARY AND CONCLUSIONS 



C. Determining 9i3 by using the day-night 
asymmetry 

As we observed earlier, the day-night asymmetry can, 
in principle, be used for determining the mixing angle 
013 if the experimental uncertainties in Am^, 0i2, and 
the day-night asymmetry itself were known to a larger 
accuracy. We now ask how small the above uncertain- 
ties must be to obtain a reasonably low uncertainty in 
^13. To estimate the uncertainty 60i3 in ^13, we use the 
pessimistic expression 



S0 



13 



9913 



I dei3 



SAn-d 



(39) 



Let us suppose that, some time in the future, we have 
determined the best-fit values for the parameters Am^ 



and 012 to be Ar 



6.9 • 10-5 eV^ and 0i2 = 33.2° 



with the uncertainties SAr 



2.5 • 10-6 eV^ and 



6012 — 0.5°. Furthermore, suppose that we have mea- 
sured the day-night asymmetry for the ES reaction to 
be An-d — 0.03 with an uncertainty SAn-d — 0.002. 
These uncertainties roughly correspond to one-tenth of 
the uncertainties of today's measurements. Next, we es- 
timate the partial derivatives of Eq. H39|l numerically and 
obtain 0013/ d An-d - -3000°, d0i3/d0i2 ~ 0.4, and 
d0iz/dAm^ ~ -2 • 106°/eV^ This gives an estimated 
error in ^13 of 50i3 ~ 12° with 6.6° from the uncer- 
tainty in An-d, 5.2° from the uncertainty in Am? and 
0.2° from the uncertainty in 0i2. (Using the uncertain- 
ties 6 An? ~ 0.8 • 10^6 eV^ and (56*12 4° suggested in 
Ref. [33, we obtain the uncertainty ^^13 ~ 10°.) Thus, 
the uncertainty in ^13 is less dependent on the uncer- 
tainty in 012 than the uncertainties in An-d and Aw?. 
We observe that the precise value of the partial deriva- 
tives depend on the point of evaluation, but that the 



We have derived analytical expressions for the day Ps 
and night Pse survival probabilities for solar neutrinos 
in the case of three neutrino flavors. The analytical result 
has been used to numerically study the qualitative effect 
of a non-zero ^13 mixing on the day-night asymmetry 
An-d at detectors. The regenerative term in the three 
flavor framework was found to be 



Pn-d = 



-cl^Ds,^^ sir? (2012) sir? (aL) 



sin^(26»i2) 



X sm 



(40) 



where K = Arr? / (2E) and the approximation is good 
for all parts of the LMA region. We have also noted that 
we regain the two flavor expression for the regenerative 
term when ^13 —> 0. An important observation is that 
the quantitative result depends linearly on the effective 
Earth electron density number used to calculate the re- 
generative term. We have argued that the value to be 
used is that of the Earth's crust. 

In the study of the day- night asymmetry at detectors, 
it is apparent that the relative effect of a non-zero ^13 
in the LMA region is increasing for increasing Arr? and 
decreasing ^12 and measured electron energy T. The de- 
pendence on 012 is also larger for smaller T and larger 
Am^. This result holds for both the elastic scattering 
and charged-current detection of neutrinos. 

We have also shown that the effects of a non-zero 6*13 on 
the isocontours of constant day-night asymmetry An-d at 
detectors are small compared to the current experimental 
uncertainty in An-d- However, should this uncertainty 
and the uncertainties in the fundamental parameters ^12 
and Am^ become much smaller in the future, the day- 
night asymmetry could be used to determine ^13 as an 



8 



alternative to future long baseline and reactor experi- 
ments. 



Acknowledgments 

We would like to thank Walter Winter for useful dis- 
cussions and Thomas Schwetz for providing the data 
of the allowed 6'12-Am^ parameter space. This work 
was supported by the Swedish Research Council (Veten- 
skapsradet), Contract Nos. 621-2001-1611, 621-2002-3577 
(T.O.), 621-2001-1978 (H.S.) and the Goran Gustafsson 
Foundation (Goran Gustafssons Stiftelse) (T.O.). 



APPENDIX A: THE CASE OF TWO NEUTRINO 
FLAVORS 

In the two flavor case, we use the parametrization 

^ : I , (Ai) 



where s = sin 6 and c — cos 9 for the leptonic mixing 
matrix. In this case, we can also write Pse as a function 
of Ps and P2e only such that 

Pn-d = PsE -PS= ^—S-iP2e " siu^ 0). (A2) 

cos(26') 

For ki and ^2, a general parametrization is 

ki = , k2 = . (A3) 

If we assume that all transitions among the matter eigen- 
states are incoherent and/or negligible in magnitude, 
then ki and fc2 are also given by 

i-Rq 

k,^ drf{r)[cos^0{r)P^^+smH{r)P^^], (A4) 

Jo 

where i € {1,2}, P^^ is the probability of a neutrino in 
the mass eigenstate |i^M,fc) at position r to exit the Sun 
in the state li^i), is the mixing angle in matter, and 
f(r) is the normalized spatial distribution function of the 
production of neutrinos. Using this, we obtain D21, as 



D2u^ / dr/(r)cos[20(r)](l-2Pju„,p), (A5) 

where fjump = ^12 — ^2^1- From this expression, it is 
easy to find that \D2iy\ < 1, and thus, < fci < 1, which 
is obviously a necessary condition. For the LMA solution 
of the solar neutrino problem, we obtain -Pjump < 10^^^"° 
by using a linear approximation of the effective potential 
at the point of resonance. This is clearly negligible. The 
value of cos(2^) is calculated as 



cos(26i) = 



Xcos(26i) - Vcc 



[Kcos{2e) - VccY+K'^ SIT? {29) 



= , (A6) 



where Vcc is the effective matter density potential. Now, 
the survival probability Ps takes the simple form 



Ps = 



1 + D2u cos(26i) 



(A7) 



For the propagation inside the Earth, we make the 
approximation that the neutrinos traverse a sphere of 
constant electron number density. This approximation 
is motivated by the fact that for current detectors, most 
neutrinos do not pass through the Earth's core. As long 
as neutrinos do not pass through the core, the only major 
non-adiabatic point of the neutrino evolution is the entry 
into the Earth's crust. As we will see, the regenerative 
term Pn-d = Pse — Ps will oscillate quickly with L and 
give an effective averaging. It follows that the exact na- 
ture of the adiabatic process inside the Earth does not 
matter as it only affects the frequency of the oscillations. 
Thus, we may use any adiabatic electron number density 
profile as long as we keep the electron number density at 
entry into the Earth's crust and at detection fixed to the 
correct values. In this case, the electron number density 
at entry into the crust and at detection are the same, 
namely the electron number density in the crust. 

Keeping the electron number density constant, also the 
effective electron neutrino potential is kept constant at 
Ve = v^GfA^b, where Gf is the Fermi couphng constant 
and Ne is the electron number density. Exponentiating 
the Hamiltonian, we obtain the time evolution operator 
and may calculate the probability P2e- The result of this 
calculation is 



P2e = sin" ( 



KVe 
4a2 



si-a^{29) sin^ (aL), 



where 



~2KVECOs{2e) + Vl 



(AS) 



(A9) 



Using the above results, we calculate Pn-d and obtain 
the result 



KVe 

= -D2.—4- sm^{2d) sm^iaL). 



4a2 



(AlO) 



For the allowed parameter space, the oscillation length 
in matter is given by Lose = l/(2a) ~ l/iiT ^ 5 • 10"* m, 
which is much shorter than the diameter of the Earth. 
Thus, we will have an effective averaging of the term 
sin (aL) in Eq. (|IlO)l . 

The expression for Pn-d contains many expected fea- 
tures. For example, when Ve = then Pn-d = 0, the 
oscillation frequency is just the one that we expect, and 
Pn-d = for D21, — reflecting the fact that this cor- 
responds to fci = fc2 = 1/2. In Eq. (|A10|1 . we observe 
that the sign of Pn-d depends on the sign of D2U, which 
depends on the point of production. If we suppose that 



^jump 1, then from Eq. (jA6|) we find that D21, is neg- 
ative if Vcc is larger than the resonance potential and 
positive if Vcc is smaller than the resonance potential. If 
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the production occurs near the resonance, then D21, will is such that Vcc is larger than the resonance potential, 
be small, since 9 ~ 45° in this case. For the allowed pa- and thus, will be negative. The same is true for D^^. 
rameter space, the region of production for neutrinos 



[1] Y. Fukuda et al. (Super-Kamiokande Collaboration), 

Phys. Rev. Lett. 81, 1562 (1998), hep-ex/9807003. 
[2] Q. R. Ahmad et al. (SNO Collaboration), Phys. Rev. 

Lett. 87, 071301 (2001), nucl-ex/0106015. 
[3] Q. R. Ahmad et al. (SNO Collaboration), Phys. Rev. 

Lett. 89, 011301 (2002), nucl-ex/0204008. 
[4] Q. R. Ahmad et al. (SNO Collaboration), Phys. Rev. 

Lett. 89, 011302 (2002), nucl-ex/0204009. 
[5] S. N. Ahmed et al. (SNO Collaboration), nucl- 

ex/0309004. 

[6] K. Eguchi et al. (KamLAND Collaboration), Phys. Rev. 

Lett. 90, 021802 (2003), hep-ex/0212021. 
[7] M. B. Smy et al. (Super-Kamiokande Collaboration), 

hep-ex/0309011. 
[8] L. Wolfenstein, Phys. Rev. D17, 2369 (1978). 
[9] E. D. Carlson, Phys. Rev. D34, 1454 (1986). 
[10] A. H. Guth, L. Randall, and M. Serna, JHEP 08, 018 

(1999), hep-ph/9903464. 
[11] A. S. Dighe, Q. Y. Liu, and A. Y. Smirnov (1999), hep- 

ph/9903329. 

[12] P. C. de Holanda and A. Y. Smirnov, hep-ph/0309299. 
[13] A. Bandyopadhyay, S. Choubey, S. Goswami, S. T. Pet- 

cov, and D. P. Roy, Phys. Lett. B583, 134 (2004), hep- 

ph/0309174. 

[14] E. Lisi and D. Montanino, Phys. Rev. D56, 1792 (1997), 

hep-ph/9702343. 
[15] M. C. Gonzalez-Garcia, C. Pefia-Garay, and A. Y. 

Smirnov, Phys. Rev. D63, 113004 ' (2001), hep- 

ph/0012313. 

[16] A. M. Dziewonski and D. L. Anderson, Phys. Earth 
Planet. Interiors 25, 297 (1981). 



[17] M. Apollonio et al. (CHOOZ Collaboration), Phys. Lett. 

B466, 415 (1999), hep-ex/9907037. 
[18] K. Hagiwara et al. (Particle Data Group), Phys. Rev. 

D66, 010001 (2002). 
[19] J. N. Bahcall, M. H. Pinsonneault, and S. Basu, Astro- 

phys. J. 555, 990 (2001), astro-ph/0010346. 
[20] M. Maltoni, T. Schwetz, M. A. Tortola, and J. W. F. 

Valle, Phys. Rev. D68, 113010 (2003), hep-ph/0309130. 
[21] Homepage of J.N. Bahcall, 

http: / / www .sns.ias.edu/'jnb/' 
[22] M. Apollonio et al. (CHOOZ Collaboration), Eur. Phys. 

J. C27, 331 (2003), hep-ex/0301017. 
[23] P. Huber, M. Lindner, and W. Winter, Nucl. Phys. 

B645, 3 (2002), hep-ph/0204352. 
[24] M. Apollonio et al., hep-ph/0210192. 
[25] P. Huber, M. Lindner, T. Schwetz, and W. Winter, Nucl. 

Phys. B665, 487 (2003), hep-ph/0303232. 
[26] Y. Itow et al., hep-ex/0106019. 

[27] A. S. Dighe, M. T. Keil, and G. G. Raffelt, JCAP 0306, 

006 (2003), hep-ph/0304150. 
[28] S. Ando, Y. H. Song, T. S. Park, H. W. Fearing, 

and K. Kubodera, Phys. Lett. B555, 49 (2003), nucl- 

th/0206001. 

[29] SNO HOWTO kit, "http://www.sno.phy.queensu.ca/sno/" 
[30] A. Bandyopadhyay, S. Choubey, and S. Goswami, Phys. 

Rev. D67, 113011 (2003), hep-ph/0302243. 
[31] For the neutral-current (NC) reaction, Vx+d — + p+n+i^x, 
the cross sections for all neutrino flavors are the same to 
leading order in the weak coupling constant. As a result, 
there will be no day-night effect in the NC reaction. 



